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Abstract 

The reactions 7p — ► K + A and 7p — ► K* + A are analyzed within perturba- 
tive QCD, allowing for diquarks as quasi-elementary constituents of baryons. 
The diquark-model parameters and the quark-diquark distribution amplitudes 
of proton and Lambda are taken from previous investigations of electromag- 
netic baryon form factors and Compton-scattering off protons. Unpolarized 
differential cross sections and polarization observables are computed for dif- 
ferent choices of the K and K* distribution amplitudes. The asymptotic form 
of the K distribution amplitude (oc £1X2) is found to provide a satisfactory 

description of the K photoproduction data. 
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I. INTRODUCTION 



The reactions 7p — > KY (Y=A,E) belong to the most elementary processes which allow 
to study strangeness production. Stimulated by the advent of a new generation of intermedi- 
ate energy electron facilities, like ELSA or CEBAF, they recently received renewed interest. 
Traditional hadronic models applied to the analysis of K photoproduction mostly make use 
of Feynman-diagram techniques The corresponding reaction mechanism is based on 

the exchange of p, A, S, K, and K*, along with a varying number of N* and Y* resonances. 
Apart from some problems with SU(3) bounds on the hadronic coupling constants (?kyn [|J) 
such models seem to work properly for photon energies up to p7 ab ^1.4 — 2.2 GeV. New data 
of higher precision and completeness (which include also spin observables) |||| are expected 
to restrict still persisting uncertainties in the meson-baryon couplings and the resonance 
parameters. 

A more fundamental treatment of photoproduction should, of course, rely on QCD, the 
dynamics of interacting quarks and gluons. A step in this direction are effective, "QCD- 
inspired" models which include already one or the other feature of QCD. A particular ex- 
ample is the chiral quark model which has been applied to K photoproduction very recently 
0] . Its elementary degrees of freedom are constituent quarks and the members of the (low- 
est lying) pseudoscalar meson octet (K, rj, and 7r). The latter are considered as Goldstone 
bosons associated with the spontaneous breaking of chiral symmetry. The Goldstone bosons 
couple directly to the (confined) quarks . With less parameters than hadronic approaches 
this model also provides reasonable results for p/ab^ GeV. Its restricted range of validity 
is caused by the use of non-relativistic transition operators and baryon wave functions. 

Direct application of QCD is (till now) restricted to kinematical situations in which the 
scattering of the hadronic constituents and their hadronization takes place on rather differ- 
ent scales. In general this means large energies and momentum transfers (p_i_). Well beyond 
the resonance region (p^ b ^> 1 GeV) exclusive photoproduction cross sections exhibit a 
characteristic angular dependence. At forward (small t) and backward (small u) angles the 
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strong variation of the differential cross section is adequately reproduced by the exchange of 
meson and baryon Regge trajectories, respectively ||. Around 8 cm = 90° (large t and u) the 
cross section flattens and shows (for fixed angles) an energy dependence typical for a hard 
interaction between the photon and the constituents inside the proton. A constituent scat- 
tering model for high-energy, large-p^ elastic and quasielastic reactions has been proposed 
in Ref. ||. The interaction mechanism of this model, namely quark interchange, may be 
also thought of as one of the simplest ways to describe photoproduction of open strangeness. 
The resulting interchange amplitude is just a convolution over light-cone wave functions of 
the interacting quarks, which have to be parameterized in an appropriate way. 

A more subtle picture, often called the "hard-scattering approach " , emerges if one tries 
to figure out the leading twist contributions to hard exclusive processes within perturbative 
QCD ||10|| . The outcome of such an analysis is a factorization formula which is also expressed 
as a convolution integral. This integral now consists of distribution amplitudes (DAs) and 
a hard scattering amplitude. The process dependent hard scattering amplitude is pertur- 
batively calculable and represents the scattering of the hadronic constituents in collinear 
approximation. The process independent DAs contain the non-perturbative bound-state 
dynamics of the hadronic constituents. DAs are, roughly speaking, valence Fock-state wave 
functions integrated over the transverse momentum. At present, the knowledge on hadron 
DAs is still rather limited. The main information is provided by QCD sum-rule techniques 
which give estimates of the lowest moments of various meson and baryon DAs |TTHl3[ . The 
few lowest moments impose some restrictions on the shape of the DAs but do not deter- 
mine them uniquely. A thorough discussion on how to construct model DAs reproducing a 



certain number of moments can be found in Ref. [14]]. One should also note that the DAs 
constrained by QCD sum rules are subject to severe criticism [p~5HT8| . Thus one is forced, at 
present, to supplement the lack of theoretical knowledge on DAs by some input coming from 
experiment. Photoproduction reactions are, in this respect, certainly very interesting. They 
exhibit a rich flavour structure and are still simple enough to allow for the computation of 
all the Feynman diagrams which enter the hard scattering amplitude. Perturbative QCD 



predictions for various photoproduction channels have been published in Ref. [ 19| . This 
paper discusses also the sensitivity of the results on the choice of the hadron DAs. Apart 
from the fact that there are objections to the numerics of this work (cf. Sect. [TV]) and it 
still needs confirmation, the predictions for the K-A channel occur to be in considerable 
disagreement with experiment. 

The present work concentrates on photoproduction of the K-A and K*-A final states. We 
consider these reactions within a particular version of the hard-scattering approach in which 
baryons are treated as quark-diquark systems. The same approach has already been applied 
successfully to other photon-induced hadronic reactions like magnetic and electric baryon 



form factors in the space- pO] and time-like region |21], real and virtual Compton scattering 



|22|| , and two- photon annihilation into proton-antiproton Further applications of the 



diquark model include the charmonium decay r] c — > pp [21] and the calculation of Landshoff 



contributions in elastic proton-proton scattering . The introduction of diquarks does not 



only simplify computations, but is rather motivated by the requirement to extend the hard- 
scattering approach from large down to intermediate momentum transfers GeV 2 ). 
This is the momentum-transfer region where experimental data are still available, but where 
still persisting non-perturbative effects prevent the pure quark hard-scattering approach to 
become fully operational. Diquarks may be considered as an effective way to cope with such 
non-perturbative effects. It is an assumption that, on an intermediate momentum-transfer 
scale, two of the three valence quarks in a baryon make up a diquark cluster. However, from 
many experimental and theoretical approaches there have been indications suggesting the 
presence of diquarks. For instance, they were introduced in baryon spectroscopy, in nuclear 
physics, in jet fragmentation and in weak interactions to explain the famous A J = 1/2 rule. 
Diquarks also provide a natural explanation of the equal slopes of meson and baryon Regge 
trajectories. For more details and for references, see f2"4]| . It is important to note that QCD 
provides some attraction between two quarks in a colour {3} state at short distances as is 
to be seen from the static reduction of the one-gluon exchange term. Also the instanton 
force seem to lead to diquark formation [^5]. Even more important for our aim, diquarks 



have also been found to play a role in inclusive hard-scattering reactions. The most obvious 
place to signal their presence is deep inelastic lepton-nucleon scattering. Indeed the higher 
twist terms, convincingly observed |26|, can be modeled as lepton-diquark elastic scatter- 
ing. Baryon production in inclusive p-p collisions also clearly reveals the need for diquarks 
scattered elastically in the hard elementary reactions P7| ]. For instance, kinematical depen- 
dencies or the excess of the proton yield over the antiproton yield find simple explanations 
in the diquark model. No other explanation of these phenomena is known as yet. 

The main ingredients of the diquark model are baryon DAs in terms of quarks and 
diquarks, the coupling of gluons and photons to diquarks, and, in order to account for the 
composite nature of diquarks, phenomenological diquark form factors. The proper choice of 
the diquark form factors guarantees the compatibility of the diquark model with the pure 
quark hard-scattering approach in the limit p± — > oo. In so far the pure quark picture of 
Brodsky-Lepage and the diquark model do not oppose each other, they are not alternatives 
but rather complements. The model parameters have been determined in Ref. [^(J by means 
of elastic electron-nucleon scattering data. The full model incorporates scalar (S) and vector 
diquarks. Vector diquarks are important for the description of spin observables which violate 
hadronic helicity conservation, i.e. quantities not explicable within the pure quark hard- 
scattering approach. The nice and simplifying feature of the two photoproduction reactions 
we are interested in is that they are not influenced by vector diquarks since only the S[ U) d] 
diquark is common to proton and A. 

The following section starts with an outline of the hard-scattering approach with (scalar) 
diquarks. It contains also a description of the Kaon, proton, and A DAs to be used in the 



sequel. Section [111] deals with the constituent kinematics, photoproduction observables, 
and the general structure of photoproduction helicity amplitudes within the diquark model. 
Predictions for photoproduction observables with a discussion of their dependence on the 
choice of the Kaon DA can be found in Sec. [TV]. Conclusions and prospects are given in 
Sec. |V|. Analytical expressions for the helicity amplitudes are tabulated in Appendix [A|. Our 
numerical method for treating propagator singularities is sketched in Appendix [B]. 



II. HARD SCATTERING WITH DIQUARKS 



Within the hard-scattering approach a helicity amplitude M^x} for the reaction 7p 
K(*)+A is (to leading order in l/p_i_) given by the the convolution integral [jlO] 



M {x} (s, t) = J dx 1 dy 1 dz 1 <p Ki * ) \z 1 , Q)4> A \yi, Q)f {X} (x lf y u z x - s, t)0 p (xi, Q) . (2.1) 

The distribution amplitudes H are probability amplitudes for finding the valence Fock state 
in the hadron H with the constituents carrying certain fractions of the momentum of their 
parent hadron and being collinear up to a factorization scale Q. In our model the valence 
Fock state of an ordinary baryon is assumed to consist of a quark and a diquark (D). We fix 
our notation in such a way that the momentum fraction appearing in the argument of 0h is 
carried by the quark - with the momentum fraction of the other constituent (either diquark 
or antiquark) it sums up to 1 (cf. Fig. [I]). In what follows we will neglect the (logarithmic) 
Q dependence of the DAs since it is of minor importance in the restricted energy range 
we will be interested in. The hard scattering amplitude T{a} is calculated perturbatively 
in collinear approximation and consists in our particular case of all possible tree diagrams 
contributing to the elementary scattering process 7uS — > ussD. A few examples of such 
diagrams are depicted in Fig. || The subscript {A} represents the set of possible photon, 
proton and A helicities. We have written the Mandelstam variables s and t with a hat to 
indicate that masses are neglected during the calculation of the hard scattering amplitude. 
They are only taken into account in flux and phase-space factors. 

If one assumes zero relative orbital angular momentum between quark and diquark and 
takes advantage of the collinear approximation (p q = x\Pb and po = ^2Pb = (1 — ^i)pb) the 
valence Fock-state wave function of a baryon B belonging to the energetically lowest lying 
octet may be written as 

* b (Pb; A) = xs «(pb, A) + /« <ft( Xl ) xv 4 {l a + —) iMPB, A) . (2.2) 



The two terms in Eq. fl2.2p represent configurations consisting of a quark and either a 
scalar or vector diquark. The pleasant feature of the covariant wave-function representation 
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Eq.( [2.2|) is that it contains, besides x\ and a (the Lorentz index of the vector-diquark 
polarization vector), only baryonic quantities (momentum p-g, helicity A, baryon mass tub)- 
For an SU(6)-like spin-flavour dependence the flavour functions % for proton and A take 
on the form (the notation should be obvious) 

Xs = uS M] , Xy = [uV {Ujd } - V2dV {U)U} ]/\/3, (2.3) 



Xt = [uS [d)S] - dS M - 2sS [U)d] ]/V6, xv = [uV {djS} - dV {u>s} ]/V2. (2.4) 

Similarly, also the q-q wave functions of pseudoscalar (PM) and vector (VM) mesons 
may be represented in a covariant way 

* PM (Ppm) = / PM PM (xx) X PM ^( A>m + m PM )75 , (2.5) 

^ VM (Pvm; A) = -/ VM VM (x 1; A) X VM ^|(^vm + m VM ) /(A) , (2.6) 

with the flavour function of the K < -*- )+ meson given by 

X KW+ =us. (2.7) 

At this point we are already in the position to recognize a considerable simplification 
in the treatment of the reaction 7p — > A as compared to arbitrary photoproduc- 
tion processes. Photoproduction of the KW -A final state can solely proceed via the S[ Ujd ] 
diquark. This is the only kind of diquark occurring in both, the proton and the A wave 
function (cf. Eqs. ( p^3|) and ( |2.4|) ). The opposite situation, namely that only the V { Uid } 
diquark becomes involved, holds for 7p — > K^*^ + S . The fact that scalar diquarks as well 
as (massless) quarks do not change their helicity when interacting with a gluon imposes 
already strong restrictions on spin observables of the K + -A and K* + - A channels. Helicity 
amplitudes which require the flip of the baryonic helicity are predicted to vanish, e.g., for the 
7P — > K + A process. On the other hand, helicity flips may take place in the 7p — > K + S° 
reaction by means of the vector diquark. In order to work out the different features of scalar 



and vector diquarks a comparison of the A and S° photoproduction channels would certainly 
be of great benefit. 

The complicated, non-perturbative bound-state dynamics is contained in the DAs H . 
These are light-cone wave functions integrated over transverse momentum (up to Q). The 
r = values of the corresponding configuration space wave functions are related to the 
constants / H . We will check the sensitivity of our photoproduction calculation on the shape 
of the K (K*) DA by choosing two qualitatively rather different forms: the one is the 
asymptotic DA 



y asy 



x) = 6x(l-x), (2.8) 



which solves the Q evolution equation for cf)(x,Q) in the limit Q — > oo (see, e.g., Ref. | 
the other one is a two-humped DA, namely 



^ z (x) = N K asy (a;)[O.O8 + 0.6(1 - 2a;) 2 + 0.25(1 - 2xf] , (2.9) 



for the K, 



(j)f{x) = A£*0 asy (x)[O.18 + 0.1(1 - 2xf + 0.41(1 - 2xf] , (2.10) 
for the longitudinally polarized K*, and 

4%*(x) = A^VasyO) [0-284 + 0.07(1 - 2x) - 0.534(1 - 2a;) 2 + 0.21(1 - 2a;) 3 + 0.267(1 - 2a;) 4 ] 

(2.11) 

for the transversially polarized K*. The DAs, Eqs. ( [2.9| )- (|2.11| ), have been proposed 
in Refs. |ll] and p9| . They reproduce the corresponding QCD sum-rule moments at 



Q 2 = 0.25 GeV 2 . It has been demonstrated quite recently that the linear x-dependence 
of pseudoscalar meson DAs at the end points x — > 0, 1 can be considered as a direct conse- 
quence of qcd g. 

The usual normalization condition, Jq dxcj) }i (x) = 1, fixes the constants N in Eqs. ( |2.9|) - 
( p.ll| ). The quantities / PM and / VM showing up in Eqs. (|2.5|) and (p.6|) are related to 
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experimentally determinable decay constants of the corresponding mesons. From the K + — > 
yU+z^t decay one infers in particular that f K+ = f^ CSLJ /2^ = 32.6 MeV. The value of f K * 
is only known indirectly via the QCD-sum-rule result f K * = 1.05f p (cf. Ref. ||11|| ). The 
experimental value of f p = /d ecay /2v^6 = 40.8 MeV obtained from the p° — > e + e~ decay 
implies f K * = 42.9 MeV. 



In previous applications of the diquark model fl20H22f a DA of the form 

<f>s{x) = N s x(l - x) 3 exp 



b 2 ( m l , m l 



X 



(2.12) 



proved to be quite appropriate for the quark-scalar diquark Fock state of octet baryons B 
The origin of the DA, Eq. fl2.12|) , is a nonrelativistic harmonic-oscillator wave function |3T| 



Therefore the masses appearing in the exponentials have to be considered as constituent 
masses (330 MeV for light quarks, 580 MeV for light diquarks, strange quarks are 150 MeV 
heavier than light quarks). The oscillator parameter b 2 = 0.248 GeV -2 is chosen in such 
a way that the full wave function gives rise to a value of 600 MeV for the mean intrinsic 
transverse momentum of quarks inside a nucleon. Note, that the DA, Eq. ( |2.12| ), exhibits a 
flavour dependence due to the masses in the exponential. The exponential in ( |2.12| ) is merely 



introduced for theoretical purposes (e. g. in order to suppress the soft end-point regions). In 
the actual data fitting the exponential plays only a minor role. Therefore, the masses and 
the oscillator parameter are not considered as free parameters but taken from the literature. 
We stress that the constituent masses do not appear in the hard scattering amplitudes. 

The dynamics of diquarks is governed by their coupling to gluons and photons. With 
respect to colour the diquark behaves like an antiquark. In order that the diquark in combi- 
nation with a colour-triplet quark gives a colourless baryon it has to be in a colour antitriplet 
state. The colour part of the quark-diquark wave function (omitted in Eq. ( |2.2j ) is therefore 
VqD° ur = (l/>/3) Ha=i S aa - The Feynman rules of electromagnetically interacting scalar di- 
quarks are just those of standard scalar electrodynamics |32|. Replacement of the electric 



charge e es by —g s t a } with g s = y/47ra s denoting the strong coupling constant and t a = X a /2 
Gell-Mann colour matrices, yields the corresponding Feynman rules for strongly interacting 
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scalar diquarks. The explicit expressions for 7-S and g-S vertices read: 



S7S : - ie e s (pi +p 2 )n , l^gS : -2ie e s g s t a g^ 



SgS : ig s t«( Pl + p 2 ), , gSgS : ig 2 s {t\ t b } 9i 



flU ■ 



(2.13) 

During the calculation of Feynman diagrams diquarks are treated as point-like particles. 
The composite nature of diquarks is taken into account by multiplying the expressions for 
the various Feynman diagrams with diquark form factors 



Qs + Q 



2 



1 n = 1 

(2.14) 

o-s n > 2 



which depend on the number (n) of gauge bosons going to the diquark. This choice of the 
form factors ensures that the scaling behaviour of the diquark model goes over into that 
of the pure quark model in the limit p± — ► 00. The factor <5 S = a s (Q 2 )/a s (Ql) (5s = 1 
for Q 2 < Qg) provides the correct powers of a s (Q 2 ) for asymptotically large Q 2 . For the 
running coupling constant a s the one-loop result a s = 127r/25 ln(Q 2 /AQ CD ) is used with 
Aqcd = 200 MeV. In addition, a s is restricted to be smaller than 0.5. The possibility of 
diquark excitation and break-up in intermediate states where diquarks can be far off-shell 
is taken into consideration by means of the strength parameter ag. 

Due to the reasons already mentioned, vector diquarks do not show up in the reactions we 
are investigating in the present paper. However, they have been dealt with in Refs. [p0| -|22 
to which we refer for further details of the diquark model. 

III. PHOTOPRODUCTION OF MESONS - 
KINEMATICS AND HELICITY AMPLITUDES 

Exclusive photoproduction of pseudoscalar mesons can, in general, be described by 4 
independent helicity amplitudes. Following the notation of Ref. [^3| we denote these ampli- 
tudes by 
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N = M 0i _ |j+li+ i , Si = M M)+1) _i , 

D = M 0>+ i i+lj _i , S 2 = M 0i+ i j+li+| . (3.1) 



N , Si, 5*2, and .D represent non-flip, single-flip, and double-flip amplitudes, respectively. 
Our helicity amplitudes are normalized in such a way that the unpolarized differential cross 
section is given by 

da 1 



dt 32vr(s - ml) 2 



(|Af + |S 1 | 2 + |,S 2 | 2 + | J D| 2 ) . (3.2) 



As we have argued above two of the four amplitudes vanish, N = D = 0, if we concentrate 
on the particular process 7p — > K + A and treat it within the diquark model. Out of the 
15 polarization observables discussed in Ref. |33| there are only 3 observables which remain 
nonzero and which differ from each other (and from da/dt) for vanishing N and D. These 
can be chosen as the photon asymmetry 

4 = ^ - ^ = 1B i 1 ™ * - ™ - > • < 3 - 3 > 

at dt dt 167r(s — mt) z 
and the two double-polarization observables 

G % = ~ Tfi ( 1 2^ + ND*) , (3.4) 
dt 167r(s — m~) J 

and 

E^r = w ( 1 2 , 2 (l^l 2 - l^il 2 + \S 2 \ 2 ~ \Df) . (3.5) 
dt 327r(s — m^y v ' 

da± (da\\) denotes the cross section for photons polarized perpendicular (parallel) to the 
reaction plane. 

Photoproduction of vector mesons may be expressed by altogether 12 linear independent 
helicity amplitudes |J33[. The diquark model leaves four amplitudes nonzero if applied to 
the formation of the K* + -A final state. In addition to Si and 5 2 one can choose, e.g., 
M + i _i +1 + i and M_ x + i +1 _i . Due to the lack of experimental data we will restrict our 
discussion of K* + -A production to the unpolarized differential cross section which is obtained 
from Eq. ( p.2[ ) by including also |M +1 _i | 2 and |M_ lj+ i )+1 _i | 2 
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The hard scattering amplitude T{a} for the elementary process 7uS[ UjS ] — > ussS[ U)S ] con- 
sists, in general, of 79 different tree diagrams. However, only 63 diagrams are encountered 
if the outgoing meson is a K + or K* + . The other 16 diagrams require the s-s pair to go into 
the produced meson. Diagrams of this type would, e.g., be important in 7p — > 0p. 

The helicity structure of the hard scattering amplitude T{\y is particularly simple for the 
K + -A and K* + -A final states. Assigning helicity labels to the hadronic constituents as in 
Fig. [I] one finds (with the S diquark helicities A2 = A6 = 0): 

A p = Ai= A3 , 

Aa = A 5 =-A 4 . (3.6) 

Thus the quark helicities are uniquely determined by the proton and A helicity, respectively. 
The additional relation (hadronic helicity conservation) 

A 3 + A 4 = A p — Aa = A K M (3.7) 

is the condition for the hard scattering amplitude T^x} and consequently the hadronic am- 
plitude M{a} to become nonzero within the diquark model. 

As depicted in Fig. |2]the hard scattering amplitude T{\y(xi, y±, Z\\ s, i) for the elementary 
process 7uS[ Uj d] — > ussSr U) d] can be decomposed into 3-, 4-, and 5-point contributions 

T{x}(xx, y u zi, s, i) = --^ e u (T ( { tf{xi, yi, z x \ s, i) + f^\x 1 , y u z x \ s, if) 

_ 77| e ud (^^(x^y^z^s,^ +fjffi(xi,y 1 ,z 1 -,s,fy , (3.8) 

depending on whether one, two, or three gauge bosons go to the diquark. The additional 
superscripts q, q, and S occurring in Eq. (|3.8| ) indicate whether the photon couples to the 
u quark, the s quark, or the S diquark, respectively. For the numerical evaluation of the 
convolution integral Eq. (|2.1| ) it is advantageous to further subdivide the various n-point 
contributions into two parts which differ by their propagator singularities: 
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{A} ^' yi ' Zl ' fl ' tJ (gl + ie)(^ + Z6') + (qi + ie)(ql + ie') ' 

^( 3 ,q), . a _ f{ 3 x}\xi,yi,zi,s,i) gfxfixuyuz^Sji) 

{A} ^' yi ' Zl ' a ' tJ (ql + ieM+ie') + (g 3 2 + <e) ( 9 ! + ie') ' 

^(4,q)/ .~t\- 4tfVi^fijjj) g{t' } q) (gi,yi,-2i;s,t) 

^ I A> I 2 ' 1 ' 2/1' ^li S ' ^7 / 2 i \/r-i2 i ' t\ / 2 i • \ ' 

^( 4 ,q)/ i\- f{\}\ x ^y^ z ^s,i) gf x f(xi,yi,zi;s,i) 

^fX\ l^l' 2/l' ^1) S ' V / 2 i ■ \ / r^2 i • A / 2 , ■ N ' 

{i (gi + ie)(Di + ie') {gi + ie) 

{A} W'V^W- (ql + ie)(g 2 +ie') + (q 2 + ie)(g 2 + ie') ' 
^( 5 ,s)/ i- ?\ f{\f\xi,yi,zi;s,t) ^g } s) (xi,yi,zi;s,t) 

r {A} (^yi^i;^*)= (Di2 + ie)(D 2 + ie0 + ■ 0.9) 

Apart from g^ 2 the q~ 2 , D~ 2 , and g~ 2 denote just those quark, diquark, and gluon propaga- 
tors which can go on-shell when integrating over xi, yi, and Z\. In order to make symmetry 
properties of the functions / and g with respect to interchange of Mandelstam variables and 
momentum fractions more obvious (cf. Appendix we have also extracted the non-singular 
gluon propagator g% 2 . Explicitly the propagator denominators read: 

<? 2 - 2/2-22-s + x 2 y 2 i + x 2 z 2 u , g\ = z 2 s + x 2 i + x 2 z 2 u , 
<?3 = 2/2^i s + x 2 y 2 i + x 2 z x ii , g\ = y x s + xtfj + x x u , 
q\ = y\z 2 s + xiy x i + x x z 2 u , g\ = y 2 z x s + y 2 i + z x u , 
q\ = y i z l s + x x y\i + x l z 1 u , D\= y x z 2 s + x 2 y x i + x 2 z 2 u , 

D 2 = y2 z \S + Xiy 2 i + x\Z\U . (3.10) 

As already indicated in Eq. Q3.9|) propagator singularities are treated by means of the usual 
ie prescription. 

Analytical expressions for the functions / and g (cf. Appendix |AT) have been derived with 
the help of "FeynArts" and "FeynCalc" [^] - two program packages written in "Math- 
ematica "which serve the automatic generation and evaluation of Feynman diagrams. Since 
"FeynArts" only contains the Feynman rules of the Standard Model it had to be extended 



to deal with S diquarks as well. The spinor techniques developed by Kleiss and Stirling |37 
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have been utilized to convert strings of gamma matrices sandwiched between spinors into 
traces which can be handled by "FeynCalc". A strong indication for the correctness of our 
results is already the agreement (apart from a detected sign error) with a previous indepen- 
dent calculation performed with "FORM" |39|] (another symbolic computer program for 
high-energy physics). Further checks of our analytical results were carried out by testing 
the U(l) gauge invariance with respect to the photon and the SU(3) gauge invariance with 
respect to the gluon. The proof of gauge invariance is facilitated by observing that not only 
the sum of all 63 tree diagrams gives a gauge invariant expression, but rather each of the 
functions / and g in Eq. ( |3.9|) is by itself gauge invariant. In addition to the gauge invariance 
tests a few diagrams were recalculated by hand. 



IV. NUMERICAL RESULTS 



Our numerical studies are performed with the set of diquark-model parameters 



f s = 73.85 MeV, Q\ = 3.22 GeV 2 , a s = 0.15 , 



(4.1) 



which has been found by fitting elastic electron-nucleon scattering data pC§ and which pro- 



vides also reasonable results in other applications of the diquark model pi| , P2| . A detailed 
explanation how the convolution integral, Eq. (|2.1| ), for the various n-point contributions 
has been treated numerically is given in Appendix |B[ At this point we only want to empha- 
size that propagator singularities have been carefully separated and integrated analytically. 
The remaining integrals could be performed by means of rather fast fixed-point Gaussian 
quadrature. 

One of the characteristic qualitative features of perturbative QCD predictions is the fixed- 
angle scaling behaviour of cross sections. Within the diquark model the 7p — ► K**) + A cross 
section behaves at large s like 



da 
~dt 



oc s 



-5 



F s (3) (- <x 2 ><y 2 >i) ~h(t/s) 



— > S 



h{t/s) . 



(4.2) 
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< %2 > and < 1/2 > denote average values of the longitudinal momentum fraction of the 
diquark in a proton or A, respectively. Equation ( (4.2|) shows that the scaling behaviour of 
the pure quark hard-scattering model [lOj is recovered in the limit s — > oo. However, at 

(3) 

finite s, where the diquark form factor F s becomes operational and diquarks appear as 
nearly elementary particles, the s~ 7 power-law is modified. Additional deviations from the 
s~ 7 decay of the cross section are due to logarithmic corrections (hidden in the functions 
.Fg and h) which have there origin in the running coupling constant a s and eventually in 
the evolution of the DAs <fi (neglected in our calculation). 

A. 7P — ► K+A 

Figure |3| shows the diquark- model predictions for s 7 da/dt along with the few existing 
large-momentum transfer data |5(J and the outcome of the pure quark hard-scattering model 
|19| (long-dashed curve). Whereas the DAs of proton and A have been kept fixed according to 
Eq. ( p.!2| ) we have varied the K + DA. The solid and the short-dashed line represent results for 
the asymptotic (Eq. ( |2.8| )) and the two- humped (Eq. (|2.9| )) K + DA, respectively, evaluated 
at Ei ah = 6 GeV. The better performance of the asymptotic DA and the overshooting of 
the asymmetric DA is in line with the conclusion drawn from the investigation of the pion- 



photon transition form factor R15IJ41] where, for the case of the pion, the CZ DA is clearly 
ruled out. There, strongly end-point concentrated DAs are also overshooting the data. Our 
findings have to be contrasted with those obtained within the pure quark-model calculation 
of photoproduction |L9 , where the asymptotic forms for both, baryon and meson DAs, give 



systematically larger results than the combination of very asymmetric DAs. However, the 



numerics of Ref. |19] must be taken with some provisio. For Compton scattering off nucleons 
it has been demonstrated |^2[ that the very crude treatment of propagator singularities 
adopted in Ref. ]19j, namely keeping ie small but finite, may lead to deviations from the 
correct result which are as large as one order of magnitude. The sensitivity of our calculation 
to the choice of the baryon DAs has been checked only with respect to their end-point 
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behaviour x — > 0,1. Neglecting the exponential factor in Eq. Q2.12|) results in a slight 



reduction of the cross section, e.g. ~ 8% at 9 cm = 90° and E^ h = 6 GeV. Deviations from 
the scaling behaviour can be estimated by comparing the dash-double-dotted and the solid 
curve, which correspond to the asymptotic K + DA at E^ h = 4 and 6 GeV, respectively. 

We have also examined the relative importance of various groups of Feynman graphs and 
found the 3-point contributions to be by far the most important. 4- and 5-point contributions 
amount to ~ 5% at 9 cm = 90° and E^ ah = 6 GeV as long as only da/dt is considered. Their 
influence decreases from larger to smaller angles. Spin observables, on the other hand, are 
much more affected by 4- and 5-point contributions. 

The three non- vanishing spin observables E, S, and G are depicted in Fig. |j. Whereas 
E measures the relative strength of the two amplitudes Si and S^, £ and G are in addition 
influenced by the phase difference of these two amplitudes. To make the interplay of the two 
amplitudes Si and S2 more obvious we have also plotted their moduli and phases in Figs. || 
and [], respectively. For both choices of the Kaon DA Si is observed to be the dominant 
amplitude in backward direction. For it remains dominant over the whole angular range. 
In contrast, £2 becomes increasingly important for if one goes from backward to forward 
direction. This behaviour is clearly reflected by the double-polarization observable E. The 
phase difference between Si and S2 varies rather moderately over the whole angular range for 
(p^sy, whereas it changes dramatically for 0^z- Unfortunately, the information on the phase 
difference is hidden in the photon asymmetry £ and the double polarization observable G 
for which the dependence on the choice of the Kaon DA is not so aggravating. 

Let us recall at this point that the occurrence of nontrivial phases in photoproduction 
amplitudes is a consequence of the fact that most of the Feynman diagrams contain internal 
gluons that can propagate on mass shell in certain kinematic regions of the momentum- 
fraction space. The treatment of the corresponding propagator singularities by means of the 
usual Feynman prescription results in an imaginary contribution to photoproduction ampli- 
tudes. One may worry about the validity of perturbation theory for a freely propagating 
gluon which is expected to be modified by long-distance effects. But fortunately photopro- 
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duction belongs to a class of exclusive reactions which does, to leading order perturbative 
QCD, not require the resummation of gluonic radiative corrections (Sudakov effects). As 
has been proved in Ref. j|3] the standard factorization formula, Eq. ( |2 . 1| ) , produces already 
an infrared finite amplitude. 

Unlike S, G, and E, which have not been measured as yet, the determination of the A 
polarization P has been attempted already {yj (for more recent efforts, cf. Ref. ||). The 
way to determine the transverse polarization of the A is to detect, in addition to the K + , 
the proton coming from the weak A — > p7r~ decay. The transverse A polarization P then 
follows from the known (P dependent) angular distribution of the weak A — > p7r~ decay. 
According to the diquark model, and also the pure quark model, P is expected to vanish 
in the hard-scattering regime (tu/s 3> rn^). The present data, however, are at too small 
s and t to allow conclusions about the validity and quality of these perturbative models. 
It would be interesting to see, whether the occurrence of sizable transverse polarizations 



at p± ~ few GeV, as observed e.g. in elastic p-p scattering [15] or inclusive production 
of hyperons in p-p collisions |5BJ, continues to the 7p — > K + A process. This would be 
an indication that, besides the perturbative mechanism, non-perturbative physics (beyond 
diquarks) is still at work. 

We have also computed differential cross sections for the reaction 7n — > K°A. For 
cos((9 cm ) > they are considerably smaller than the corresponding 7p — > K + A cross 
sections. The amount of suppression depends on the choice of the K° DA. For the asymptotic 
DA the suppression factor is ~ 10 in the whole forward region, whereas it increases for the 
two-humped DA from 2 to m 10 when cos(# cm ) is varied from to 0.8 (E^ = 6 GeV). In 
view of the plans at CEBAF to study 7n — > K°A by means of a deuteron target || this is 
certainly an interesting observation which could be helpful to pin down the uncertainties of 
the Kaon DA. 
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B. 7P — ► K*+A 

The diquark-model results for (da/dt)^ p ^ K ,+ K are plotted in Fig. [/]. Again curves are 
shown for two choices of the K* + DA with p and A DA kept fixed (E^ = 6 GeV). The 
results resemble those for photoproduction of K + mesons. Cross sections for the asymmetric 
DA, Eqs. (P7T0) and ( j2.11| ), are nearly one order of magnitude larger than for the asymptotic 
DA, Eq. (|2.8| ). If asy is taken for both, K + and K* + , the photoproduction cross section for 
the K* + vector meson is found to be by a factor of 1.8 - 3.6 (depending on the scattering 
angle) larger than that for the pseudoscalar K + meson. An increase by a factor 1.73 is 
due to the different K + and K* + decay constants f K and f K *. The remaining difference 
is caused by the contribution of transversially polarized K* + mesons, which increases from 
small to large scattering angles. The situation would be quite similar if we had taken <pQ Z 
for both, K + and K* + . However, the enhancement of the K* + cross section is completely 
compensated, if is replaced by (fyf and 0^* (cf. Eqs. ( |2 . 1 0|) and ( |2.11|) ) when going from 
K + to K* + photoproduction. 

Till now large momentum transfer data for photoproduction of vector mesons are only 
available for the reaction 7p — > (p + ^)p [0. At 9 cm = 90° and ~ 6 GeV the 
experimental value of the cross section ratio (da /dt) Jp ^( po+OJ ) p / (da / dt)^^^ is ~ 2. This 
means in particular that (da/dt) 1P ^ Popi <2(da /dt) 7P _ )7ro p. The difference in the po and no 
decay constants (f^ ecay ~ l-5/J ecay ), however, already implies an enhancement of the p 
photoproduction cross section as compared to the ir Q one by a factor of ~ 2.3 which is further 
magnified by contributions of transversially polarized ps. The only way to compensate part 
of this enhancement is to assume that (as above in the case of K + and K* + ) the DAs of 7r and 
p differ from each other. Experimental data on photoproduction of pseudoscalar and vector 
mesons with the same flavour content could thus be very useful to work out differences in 
the corresponding DAs. 
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V. SUMMARY AND CONCLUSIONS 



We have investigated photoproduction of K-A and K*-A final states in the few-GeV 
momentum-transfer region. Our analysis is based on perturbative QCD supplemented by 
the assumption that baryons can be treated as quark-diquark systems. The present cal- 
culation continues previous work on photon-induced hadronic reactions PD|-f2^] performed 
within the same approach. By modeling quark-quark correlations inside a baryon as quasi- 
elementary particles - scalar and vector diquarks - we account for some non-perturbative 
effects. In this way we are able to extend the range of applicability of the pure quark hard- 
scattering approach from large down to moderately large momentum transfers. The fact 
that the photoproduction channels we are interested in contain a A in the final state entails 
a considerable reduction in computational effort. In contrast to arbitrary photoproduction 
reactions only scalar diquarks must be taken into consideration. This has the consequence 
that helicity amplitudes and hence spin-observables violating hadronic helicity conservation 
(cf. Eq. ( p.7|) ), e.g. the A polarization, are predicted to vanish. 

Our numerical studies have been performed with the diquark-model parameters and the 
quark-diquark DAs proposed in Ref. p0| . We have paid special attention to the correct 
and numerically robust treatment of propagator singularities (cf. Appendix [FJ). Reasonable 
agreement with the few existing 7p — > K + A data is achieved already with the asymptotic 
form for the K + DA. On the other hand, the end-point concentrated K + DA, based on QCD 



sum rules [0]], seems to perform less well. The corresponding curve lies far beyond the data. 
The difference between these two Kaon DAs is also clearly visible in the three non-vanishing 
polarization observables, i.e. the photon asymmetry S and the two double-polarization 
observables G and E. It is most pronounced in the observable E. Another quantity which 
we found to be very sensitive on the choice of the Kaon DA is the angular dependence of 
the cross-section ratio (da/dt)^ n ^ K o\/(da/dt)^ p ^K+A- For the photoproduction of the K* 
vector meson there are no data to compare with. We have again tested the asymptotic 
DA and a K* DA which obeys QCD sum-rule constraints on the lowest moments [29||. The 
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differences in the results for the two DAs are quite similar to those for photoproduction of 
the pseudoscalar K + meson. When going from K to K* photoproduction the cross section 
becomes larger due to the different K and K* decay constants and the additional contribu- 
tions from transversially polarized K*s. This increase of the cross section, however, is partly 
compensated if different DAs for K and K* are used. 

With regard to future experiments we consider, apart from more and better large p± cross 
section data, the polarization measurement of the recoiling A as one of the most urgent tasks. 
A large polarization indicates that the perturbative QCD regime has not been entered yet. 
In the perturbative QCD regime the Kaon DAs could be restricted by means of quantities, 
like the photon asymmetry or the cross section ratios (da/dt) jn ^ K QA/(da/dt)^ p ^ K + A and 
(da/dt) 7p ^ K *+ A /(<ia/<it) 7 p^K+A, which are very sensitive to the choice of the DAs. With a 
maximal photon laboratory energy of (at present) 4 GeV CEBAF |J touches at best the 
border of the hard-scattering domain. More decisive data could be expected from a future 
electron facility like ELFE ||7| which is designed to explore the energy range up to 15 GeV 
(or even higher) with a continuous high intensity electron beam. 



In this appendix we quote analytical expressions for those hard amplitudes which describe 



speaking, these are just the scalar diquark contributions to photoproduction of pseudoscalar 
mesons. According to Eq. Q3.9|) the various n-point contributions to these amplitudes can be 
decomposed into gauge invariant functions / and g. The functions / and g which determine 
the hadronic (helicity conserving) amplitude Si read 



APPENDIX A: ANALYTICAL EXPRESSIONS 



FOR THE HARD AMPLITUDES 



the process 7uS(ud) 



ussS(ud) with the ss-pair being in a spin-zero state. More generally 
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V3' 



Cp "73' 



(A2) 



A^' n) = 128Tr 2 y/Traa 2 s (iu/s)\/-tF^ n \-x 2 y 2 t), 



(A3) 



where a denotes the fine-structure constant. 

The functions / and g contributing to the hadronic amplitude S 2 are obtained from those 
entering Si by interchange of the Mandelstam variables s <-> u and the momentum fractions 
x\ <-> yi and zi <-> z 2 . One finds in particular 
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(A4) 

The hard amplitudes for photoproduction of longitudinally polarized K*-mesons are (up 
to a sign change in T _i +1 _i) also determined by these expressions. For photoproduction 
of transversially polarized K*-mesons the functions / and g are of similar length and shape 
and can be obtained from the authors on request. 



APPENDIX B: NUMERICAL TREATMENT OF PROPAGATOR POLES 

The numerical difficulties in performing the convolution integral Eq. (|2.1|) are mainly 
caused by the occurrence of propagator singularities in the range of integration which give 
rise to a principal value integral 

In what follows, the four cases to be distinguished will be discussed separately. 
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1. No propagator on shell 



The only contribution to the convolution integral Eq. (|2.1| ) exhibiting no propagator 
singularity is the one corresponding to gr^P (cf. Eq. ( |3.9| )). In these circumstances the 
convolution integral is easily performed by means of 3-dimensional Gaussian quadrature. 



2. One propagator on shell 

If only one of the propagators goes on shell within the integration region - this happens 
for gf\f\ ^{a}^ ' /{a'^j anc ^ 9{x} ~ the corresponding integrals over x% have the general 
structure: 

l(\y uZl) =J q dxi __ . ( B2 ) 

In order to simplify notations we have neglected helicity labels and the dependence on the 
Mandelstam variables s and i. Furthermore, the distribution amplitudes 4> p (xi), <f>\(yi), and 
4>\^{zi) have been absorbed into the function h(xi,yi, z\). The integral 1^ may be rewritten 
to give 

T (k) ( \ f 1 , K x ^v^ z i)-K x x\y^ z i) 

I () (y 1 ,z 1 )= dxi 



o 



k 2 

1 dxi 



dk 2 



dxi 



-v 



(B3) 



where (yi, zi) represents the zero of k 2 (considered as a function of X\). The 

first integral in Eq. ( |B3| ) is now again tractable by simple Gaussian quadrature, whereas the 
principal-value integral can be done analytically. Since Jw is a regular function of y\ and 
Z\ Gaussian integration can also be applied to these variables. Analytical expressions for 
propagator-pole positions, principal value integrals and Xi-derivatives of propagator denom- 
inators are listed in Table |. 
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3. Two propagators on shell — propagator poles not coinciding 

If two propagators fej" 2 and A>T 2 go on shell one can proceed similarly as in the one- 
pole case, provided the zeroes Xi = x[ kl \yi, Z\) and x^ = x^ 2 \yi,Zi) of k\ and k\ 
(considered as functions of x\) do not coincide for fixed y\ and zi, < yi, Z\ < 1 arbitrary. 
This is guaranteed for the Feynman diagrams contributing to fj^} \ /{a}> /{A} ' anc ^ #{!'}■ 
The xi-integrals to be considered have the general form 

If Xi xj[ a partial fractioning yields 

^.*) = ^rbp{(S" V,( »''*' ) -(S" ll<W( "'"* i) } ■ (B5) 

i.e., two terms which again can be treated according to Eq. (|B3| ). 

4. Two propagators on shell — propagator poles coinciding 

This is the worst case and shows up in connection with the functions gf\}\ /{A} > an d 
/na . The general structure of the Xi-integrals is again that of Eq. (|B4|) . However, now it 
happens that the two propagator singularities x^ and x^ (which still depend on y\ and 
zx) coincide for a certain value of y\ (0 < y\ < 1), < Z\ < 1 fixed. We denote this value by 
y{kiM) _ yY uk ^(zi) - it still depends on z\. We note that the partial fractioning can still be 
performed for arbitrary values of y\ and Z\ as long as e and e' are kept finite. By carefully 
taking the limit e — > in the terms containing k\ and e' — > in the terms containing fc| and 
using Eq. ( p3|) one ends up with 
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In Eq. ( p6|) e stands for either e or e' . Closer inspection of (x% 



x 



(considered 



as function of y{) reveals that the zero y^ l,k2 ' is quadratic. On the other hand one also 
finds that y^ 1 '^ is a single zero of both h(x[ kl \yi, z\) and h(x[ k2 \ y i: Zi). This immediately 
implies that the real part of J^ fcl ' fc2 - ) (y 1 , z\) is a regular function of y\ and also z x so that 
Gaussian quadrature is again applicable to the corresponding integrations. Taking further 
into account that (dk\/ dx\)~ l and (dk^/dxi) -1 have different signs the imaginary part of 
I^ kl,k2 \yi, Zi) can be written as 

-1 
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(kl) 



dxi 

(yx - y[ klM) ) 
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(B7) 



where h — (yx — y[ u )h. Integrating ^/( fcl ' fe2 ) with respect to y\ and letting e — > gives 

dy 1 ^' ha \y 1 ,z 1 ) = p f dy^l^iy^zx) , (B8) 

Jo 

i.e., only the principle- value part of the integration survives. The principle- value integral in 
Eq. ( p8|) can be treated analogous to the principle- value integrals in x\ (cf. Eq. (|B3|)). 

Proceeding along the steps outlined in this appendix, i.e. carefully separating the singu- 
lar contributions, exploiting delta functions, rewriting principal-value integrals as ordinary 
integrals plus analytically solvable principle- value integrals, it is finally possible to do all the 
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numerical integrations by means of fixed-point Gaussian quadrature. For our purposes an x- 
y-z grid of 20 x 20 x 24 turned out to be sufficient. Taking instead a 32 x 32 x 48 grid changes 
the results by less than 0.2%. The numerical calculations were performed on a DEC7000-610 
APLPHA workstation. For the larger grid size the calculation of da/dt^p — > K + A) took 
about 1 second per energy point and angle. 
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FIGURES 

FIG. 1. Constituent kinematics for 7p — > K^+A. 

FIG. 2. A few representative examples of Feynman diagrams contributing to the elementary 
process 7uS [U]d] — * ussS [Ujd] . 

FIG. 3. Differential cross section for 7p — > K + A scaled by s 7 vs. cos(# C m)- Solid 
(dash-double-dotted) line: diquark- model result at = 6GeV (4GeV), proton and Lambda DAs 
chosen according to Eq. fl2.12Q , Kaon DA according to Eq. ( |2.8D (asymptotic DA); short-dashed line: 



diquark-model result at = 6GeV, proton and Lambda DAs chosen according to Eq. ( 2.12 ) 



Kaon DA according to Eq. ( [2.S| ) (Chernyak-Zhitnitsky DA [11Q; long-dashed line: quark-model 



result [19 for the asymmetric proton and Lambda DAs of Ref. |L^] and the Kaon DA Eq. (| 



Experimental data are taken from Ref. [40]. 



FIG. 4. Diquark model predictions for the non-vanishing 7p — > K + A polarization observables. 
Full (short-dashed) line: same as in Fig. |3|. 

FIG. 5. Differential cross section for 7p — > K + A scaled by s 7 vs. cos(# cm ) at p^ ab = 6 GeV - 
contributions of the helicity amplitudes Si and S2, respectively. Solid (long-dashed) line: contri- 
bution of Si (S2) for proton and Lambda DAs chosen according to Eq. (|2.12j ), Kaon DA according 
to Eq. ( |2.8|) ; short-dashed (dash-dotted) line: contribution of Si (S2) for proton and Lambda DAs 
chosen according to Eq. ( |2.12[ ), Kaon DA according to Eq. ( |2.9| ). 

FIG. 6. Phases of the helicity amplitudes Si and S2, respectively. Lines as in Fig. 0. 
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FIG. 7. Differential cross section for 7p — ► K* + A scaled by s 7 vs. cos(6* cm ) at pZ b = 6 GeV. 
Solid line: diquark- model prediction, proton and Lambda DAs chosen according to Eq. ( p. 12 ), K* 
DA according to Eq. ( |2.8| ) (asymptotic DA); short-dashed line: diquark-model prediction, proton 
and Lambda DAs chosen according to Eq. ( |2.12 ), K* DA according to Eqs. ( |2.10| ) and ( 2. lip (taken 
from Ref. pi). 
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TABLES 

TABLE I. Propagator-pole positions, xi-derivatives, and principal value integrals for the sin- 



gular propagators occurring in Eq. (3.9). 
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